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The oscillation and stability of differentially rotating 
spherical shells: The normal-mode problem 
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ABSTRACT 

An understanding of the dynamics of differentially rotating systems is key to many 
areas of astrophysics. We investigate the oscillations of a simple system exhibiting 
differential rotation, and discuss issues concerning the role of corotation points and 
the emergence of dynamical instabilities. This problem is of particular relevance to 
the emission of gravitational waves from oscillating neutron stars, which are expected 
to possess significant differential rotation immediately after birth or binary merger. 
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1 BACKGROUND AND MOTIVATION 

This paper describes a dynamical problem that is rele- 
vant to many areas of astrophysics. We consider the os- 
cillations of a simple system exhibiting differential rota- 
tion, and discuss issues concerning the role of corotation 
points and the emergence of dynamical instabilities. This 
is known to to be a technically challenging problem (anal- 
ogous to that of Couette flow in standard fluid dynamics) 
that remains to be understood in detail. In addition to be- 
ing mathematically interesting, the problem is of potentially 
great astrophysical relevance since the oscillations of com- 
pact objects such as neutron stars or strange stars are a 
promising source of detectable gravitational waves. Of par- 
ticular interest are oscillations that can grow via either a 
dynamical instability or a secular instability to the emis- 
sion of gravitational waves, such as the Chan drasckhar- 
Friedm an Schutz (CFS) instability dFriedman fc Schut J 
If978allrl . The most promising instability scenarios con- 
cern hot newly born neutron stars, which one would ex- 
pect to possess some differential rotation. This expectation 
is borne out by the latest studies of rotational core col- 
lapse, which indicate that the remnant will indeed emerge 
from the collapse rotating differentially llZwerger fc Mulled 
Il997l: iDimmelmeier. Font fc Mulled 120021). Subsequent ac- 
cretio n of supernova remnant material CWa^^^^AnderesorJ 
2002), or material from a companion iFuiimotol 11993). 
may also drive differential rotation. Differential rotation 
may also be maintained by the oscillations themselves. 
Recent work on r-modes in neutron stars suggests that 
the modes can drive a uniformly rota ting star into diffcr- 
ential rotation via non-linear effects jRezzolla et all l2000t 
[Levin fc Ushomirskvl 1200 ll : iLindblom. Tohline fc Vallisneril 
200f). A differentially rotating massive neutron star 



may also be generated following a binary neutron star 
merger. The maximum mass of a differentially rotating 
merger remnant can exceed significantly the maximum 
mass of a uniformly rotating rem nant, remaining stable 
over many dynamical timescales JShibata fc Urvul [2000; 
IBaumgarte. Shapiro fc Shibatall2000il . 

It is clear that differential rotation is an important 
factor in modelling these systems. Nevertheless, the vast 
majority of secular instability studies to date have fo- 
cused on unifor mly rotating neutron stars (exceptions in - 
clud e studies by|ljnamura_et_aj [ ill 995l) . lYoshida et all J2002t . 
and iKarino & Eriguchin2002T) 'l. The situation is some- 
what different for the dynamical bar mode insta bility 
(see for example Shibata. Baumgarte fc Shapirol <|2000T) and 
iNew fc Shapirol |200J| '). There have been a number of large 
scale numerical simulations aimed at understanding the bar 
mode instability. These studies provide significant insights 
into the dynamical stability of differentially rotating sys- 
tems. Yet it is clear that a better theoretical understanding 
of this problem is desirable. This need is well illustrated 
by the recent discovery of dynamical instabilities that be- 
come active at very low values of T/\W\ (the ratio of ro- 
t ational kinetic energy to gravitational potential ene rgy) 
JCentrella et all200lt IShibata. Karino fc Eriguchill2002l) . 

There are two issues associated with differential rota- 
tion that make this problem particularly interesting. Firstly, 
differential rotation may introduce new dynamical instabil- 
ities jBalbinskillf98llLuvtenlll990ah . Secondly, the dynami- 
cal equations are complicated by the presence of corotation 
points. These are points where the pattern speed of a mode 
matches the local angular velocity, and at these points the 
governing equations are formally singular. This gives rise to 
a continuous spectrum and could in principle permit other 
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corotating solutions. These two important issues have how- 
ever attracted little attention in the literature. In particular, 
it is not at all clear what the dynamical role of the contin- 

uou s spectrum may be. 

iKarino. Yoshida fc Eriguchil J200ll) studied r-mode os- 
cillations in rapidly and differentially rotating Newtonian 
polytropes for two particular rotation laws. They found 
that the r-mode frequencies shifted towards the corota- 
tion region as the amount of differential rotation increased 
and at high differential rotation the codes failed, ap- 
parently due to the presence of corotation points. Solu- 
tions with corotation points were not, however, investi- 
gated. Two subsequent studies of a differentially rotat- 
ing shell for the same two rotation laws also focused 
on modes without corotation points | Rezzoll a fc Yoshidal 



l200ll: lAbramowicz. Rezzolla fc Yoshidail2002l) . The authors 
claimed that for these rotation laws the lowest frequency 
r-modes did not develop corotation points, even at high dif- 
ferential rotation. 

For papers that deal more fully with dynamical in- 
stabilities and corotation points we must go back a num- 
ber of years (or consider the related problem of differen- 
tially rotating disks for which there is a vast literature). In 
the 1980's Balbinski carried out a number of studies on z- 
independent modes of differentially rotating cylinders (for 
cylin drical coord in ates ro , <f), z) that addressed both issues 
jBalbinskilll984allbl Il985h . The first paper reported an ana- 
lytic solution for the continuous spectrum for one particular 
rotation law. Balbinski then solved the initial value problem 
for simple initial data, and found that the physical pertur- 
bation associated with the continuous spectrum died away 
on a dynamical timescale. The second paper examined the 
secular stability of the continuous spectrum to the emis- 
sion of gravitational waves. The final paper investigated the 
behaviour of modes that had no corotation points at low dif- 
ferential rotation. As the degree of differential rotation in- 
creased, some of the modes were seen to enter corotation and 
in due course merge to give rise to a dynamical instability. 
In the early 1990's, Luyten ge neralised the wo rk in lBalb inski 
( 1985) to other rotation laws iLuvtenl |l990a). In two subse- 
quent papers he solved for the same 2-independent modes 
in sel f-gravitating homoge neous masses rather than on cylin- 
ders jLuvtenll990blll99l|) . Although Luyten's method could 
not track modes with corotation points, modes appeared to 
enter and leave the corotation region as the degree of dif- 
ferential rotation increased. Luyten too reported finding an 
instability at high differential rotation. 

There is clearly a need to analyse modes with a z- 
dependence such as the toroidal r-modes. As a first step to- 
wards this goal, we present a detailed analysis of the modes 
of a differentially rotating spherical shell. We hope to im- 
prove our understanding of the technical challenge posed by 
the presence of corotation points, and develop techniques 
that may later be generalised to the full problem of dif- 
ferentially rotating stars. By investigating dynamical insta- 
bilities we may also be able to rule out certain rotation 
laws, or set a threshold for the maximum degree of dif- 
ferential rotation that could be sustained. We select four 
rotation laws from the literature in order to compare our 
results to previous work, and one new law. We note how- 
ever that in all but one case these laws were selected for 
mathematical convenience rather than physical reasons. One 



would clearly prefer to use a rotation law for a young neu- 
tron star that is grounded in physi cs, emerging for ex ample 
from simulations of core collapse llZwerger fc Mulled IT991 : 
iDimmelmeier. Font fc Mullerll2002l) . We hope to do this in 
future work. 



2 FORMULATING THE PROBLEM 

We consider an incompressible fluid confined to move in a 
thin spherical shell. We take the shell radius to be R. Its 
thickness is taken to be so small that we do not have to ac- 
count for dynamics in the radial direction. On these assump- 
tions, and taking the perturbation to depend on exp(irrup), 
the ^-component of the Euler equation simplifies to 

1 



dtSvg + imSlSve — 2Q, cos 05v^ 
and the (^-component to 
dtSv v + imQSvp + QSvg 



pR 



dgSP, 



x8P, 



pR sin 9 

where the equilibrium vorticity is 
Q = 2Qcos9 + sm6dgD.. 

The continuity equation for an incompressible fluid is 
dg(sin68vg) + im5v v = 0. 



(1) 



(2) 



(3) 



(4) 



This tells us that only toroidal modes are permitted in this 
problem. Defining the stream function (in the standard way) 
as 

U — Svg sin 0, (5) 

we can combine equations (0 and (|5J to give the vorticity 
equation 



sin o 

where 

2 fl 2 . cos0 m 2 
Vg = dg + -^^de r^rr 



(6) 



(7) 



sin f sin t) 
is the Laplacian on the unit sphere. Changing variables to 
x = cos 9 the vorticity equation becomes, 

2xd x 



(idt—mO,) 



(I - x 2 )dl 



U-md x QU = 0.(8) 



1-x 2 

In terms of x the equilibrium vorticity is 

tl — 2xfl — (1 — x 2 )d x Q.. (9) 

Equation JHJ is singular at the poles (x = ±1). Performing 
a Frobenius expansion around these points, 



U 



(I ±x) a ^a k {l ±x) 



(10) 



we find that the regular solution has a = m/2. Factoring 
out this behaviour and defining a new variable u by 

U = u{l - x 2 ) m/2 , (If) 

we find that u obeys the equation 

(mfi — idt) [(I — x 2 )d 2 — 2(m + l)xd x — m(m + I)] u 

+md x Q,u = 0. (12) 
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We will use both equations JSJ and 1121 in this paper. 



3 UNIFORM ROTATION: THE R-MODES 

Before we proceed to discuss differential rotation it is use- 
ful to recall the standard result for uniform rotation. After 
all, we are interested in understanding how the mode spec- 
trum changes as differential rotation is introduced. The os- 
cillations of a uniformly rotating incompr essible shel l have 
been studied by a number of authors (|Haurwitj |1940|; 
IStewartson &: Rickardl Il969t iLevin fc Ushomirskvl 12001). 
and are sometimes referred to as Rossby-Haurwitz waves. 
Rossby-Haurwitz waves are however the limit on the shell of 
the r-modes of stellar perturbation theory, and we will use 
the latter term. For uniform rotation we have 



d e Q 



-2 sin en. 



(13) 



Assuming a mode solution with time dependence oc 
exp(— iat) equation © becomes 



(a - mQ)V 2 e U - ImQU = 0. 
Assuming that a 7^ mQ we have 

v2 = _amn K 

a — mil 

Expanding U in spherical harmonics 



then by Legendre's equation we must have 

i(i + i)cpYr. 



I 



Thus we must have 



i(i + i)crYr ' 2 "' n .. Yc y, 



a — mQ ' 

l 



(14) 



(15) 



(16) 



(17) 



(18) 



from which we deduce that the mode solutions are such that 
(i) only a single C; m is non-vanishing and (ii) the correspond- 
ing frequency is 

2mQ 



a = mQ 



(19) 



1(1 + 1) 

These are the familiar r-modes. 

There is, however, another possible solution, the 
geostrophic solution. Consider what happens if we set a = 
mQ. It is clear from equation 1141 that the only possible 
solution in this case is U = 0, which corresponds to a zero 
velocity perturbation. Consider however what this means for 
the displacement £. In the case of uniform rotation, 

Sv = i(mQ - cr)£e, (20) 

and 



8Vip — i(m£l — er)£ ¥ 



(21) 



It is clear that, if a = mQ., a vanishing velocity perturbation 
is compatible with a non-zero displacement. This solution, 
trivial in the case of uniform rotation, must be considered 
when the rotation is differential as it is this solution that 
gives rise to the continuous spectrum. Note the similarity 
with the standard r-modes, which are a trivial solution of 
the perturbation equations for a non-rotating star. 



4 THE NORMAL MODE PROBLEM 

Let us now consider the case of a differentially rotating 
shell. Assuming a mode solution with time dependence oc 
exp(— iat) equation I12H becomes 

(mQ — a) [(1 — x 2 )u" — 2(m + l)xu — m(m + l)u\ 

+mQ'u = 0, (22) 

where the primes indicate derivatives with respect to x. At 
the poles, we require that the solution be regular. In addi- 
tion, symmetry allows us to require that u be either even or 
odd under reflection through the equator. This implies that 
we can set boundary conditions on the equator that either 
u or u' be zero there. These boundary conditions allow us 
to solve the equation on the range < x < 1 rather than on 
the full range — 1 < x < 1. 

If a 7^ mQ(x) for all points on the shell we have a 
regular eigenvalue problem. If, however, a — mQ(x) at a 
point x = x c , then x c is a corotation point and equation 
I22t is formally singular. Corotation points are points where 
the pattern speed of a perturbation is equal to the local 
angular velocity. Note that the term corotation point is a 
little misleading ELS X — O&c 

really represents a line of latitude 
on the shell. Equation 1221 has three classes of solution: 

• Discrete real-frequency modes outside corotation, i.e. 
for which a 7^ mQ(x) for \x\ < 1. 

• Discrete complex frequency modes for which a = oy + 
iai. Dynamically unstable (growing) modes correspond to 

CTi > 0. 

• Corotating solutions for which a = mQ(x) at a corota- 
tion point x c , where \x c \ < 1. 

The first two kinds of solutions are easy to determine numer- 
ically for any given rotation law. In contrast, the solutions 
with corotation points present a real challenge. Nonetheless, 
these solutions are of great interest, not least because of their 
association with the continuous spectrum. For this reason 
much of the following discussion will concern solutions of the 
third kind. In the following sections we shall refer to solu- 
tions with corotation points as being "in corotation" . Modes 
without corotation points that subsequently develop them 
as the amount of differential rotation is increased will be 
referred to as "crossing the corotation boundary" . It should 
be noted, however, that solutions in corotation are in fact 
only corotating at one particular latitude. 



5 A NECESSARY STABILITY CRITERION 

Before discussing detailed solutions for particular rotation 
laws it is useful to check whether unstable mode solutions are 
likely to exist. Such solutions would typically correspond to 
dynamical instabilities, and hence they are of great physical 
significance. 

For inviscid parallel shear flow it is possible to derive a 
necessary condition for stability that depends on the form 
of the rotation l aw, the Rayleigh infle xion point criterion 
(see for example iDrazin fc Reidl fl982'l'l. Here we derive a 
similar condition for modes of a differentially rotating shell. 
We start with equation ||SJ, assuming a mode solution with 
complex a, 
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(mQ — a) 



-\{l-x*)U') 



m 2 U 



+ mQ'U = 0. (23) 



dx LV ' ' 1 — x 2 

Multiply by the complex conjugate U* and integrate to get 



[{l-x 2 )U'U'}\- [\l-x 2 )\U'\ 2 dx 
Jo 



m I ^rdx + m 

1 — x z 



1 tr|t/ 2 | 

mQ — a 



dx = 0. 



(24) 



The first term is clearly zero due to the boundary conditions. 
The imaginary part of equation 12411 is 

f 1 Q'\U\ 2 

ma > / i o t 9 dx = °- ( 25 ) 
Jo \mQ-a\ 2 

For dynamically unstable modes cr, 7^ and equation 12511 
can only hold if Q! changes sign at some point on the inter- 
val. This is a necessary (although not sufficient) condition 
for instability. 



6 COROTATING SOLUTIONS 

We begin our study of corotating solutions by considering 
the case when there is a corotation point in the domain 
< x c < I (the symmetry of the problem means that there 
will be a second corotation point in the lower hemisphere, 
at x = —x c ). We will consider the special cases of x c = 1, 
(corotation point at the poles or the equator) later. Employ- 
ing the method of Frobenius to investigate the solutions in 
the vicinity of the corotation point, we let 



u = a n (a 



\n+7 



(26) 



For all of the rotation laws examined in this paper mQ — a 
contains only a single factor of x — x c , also Q is an analytic 
function of x. Assume also (for now) that Q' does not have 
a zero at x = x c . Substituting into equation 1221 to obtain 
the indicial equation, we find that 7 = 0, 1, so we will have 
one regular and one singular solution. The general solution 
is 



U — U re g T Using , 

where 



= a n (x — x c ) n+1 , 



(27) 



(28) 



and 



Using = [a n (x - x c )" +1 ln(x - x c ) + c„ (x - I c )*] . (29) 

Note that although u S i ng is finite and continuous at x c its 
derivatives are singular. 



6.1 Possible regular solutions 

Let us begin by asking whether the problem admits any 
solutions that are regular at corotation (i.e. contain u reg 
only). First of all note that if u is regular in corotation then 



so is U . We now apply the method that was used to derive 
a necessary condition for instability in section 

We begin with equation \1'M . multiply by U* and as 
before integrate. This time, however, the frequency is taken 
to be real and we are going to treat the two domains < 
x < x c and x c < x < 1 separately. On the first domain 



[(1 



m 



)U'U* 



It/ 2 



(l-x 2 )\U'\ 2 dx 



- x 



7 da; + m 



n'\u 2 

mQ. — 



-dx = 0. 



(30) 



'0 jo 
The first term vanishes due to boundary conditions and we 
are left with 

~ 2 \U 2 \ 



(l-x 2 )\U'\ 2 + 



rn 



&\U 2 
mQ. — 



dx + 



-dx 



It is clear that we must have 
Q'\U 2 ' 



mQ 



-dx > 0. 



(31) 



(32) 



(Note that U/(mQ — a) is well behaved at x c if u on u reg .) 
Treating the second domain in the same way we can derive 
a second condition, which is that we must have 

' 1 dx > 0. (33) 
mQ — a 

Now mQ— a clearly changes sign at the corotation point. 
So unless Q' also changes sign somewhere on the shell it 
will be impossible to satisfy both conditions. A rotation law 
will therefore only admit purely regular solutions if Q' = 
somewhere on the shell. 

There is one other class of regular solution that might be 
encountered. The governing equation 1221 will remain non- 
singular if, for a particular frequency, the zero of mQ — a 
coincides with that of Q' . This will in fact always be the 
case for some frequency in the continuous spectrum range 
if Q' has a zero on the shell. Unlike the regular Frobenius 
solution u reg , however, solutions in this class need not have 
a zero in the eigenfunction at the corotation point. 

6.2 General (singular) solutions 

The coefficients of both the regular and singular solutions 
(equations 1281 and 1291 ) are found by substituting the so- 
lutions into equation 1221 . We find Co to be related to ao, 

Q'(x c ) 



-(x c - l)o . 



(34) 



Q'(x c ) 

Hence ao and ci (equation I29|l are the two arbitrary con- 
stants in the general power series solution. 

The Frobenius expansion gives us an approximation to 
the general solution close to the corotation point. This solu- 
tion is useful in a numerical search for solutions. However, 
the logarithm term causes a problem. In particular, it gives 
rise to a branch point at x c and we would in principle need 
to integrate around this point in the complex plane. An al- 
ternative is to seek an equivalent well-behaved solution for 
x < x c and ma tch the two solu tions at the corotation point. 

Following Bal binsk] jl985h we make a change of variable 



Differentially rotating spherical shells: Normal modes 5 



2 = — . (35) 

X 

The corotation point is now at z = 1. Equation 1221 becomes 



(mil — a) 



1 2 2 

[z -x c 



\z 2 u" + 2zu) + 2(m + l)zu 



m(m + l)u 



+ muQ.' = 0. (36) 



As before, we look for a power series solution in the vicinity 
of the corotation point, letting 



u = b n (z — l) n 



+7 



(37) 



Once again the indicial equation yields 7 = 0, 1, and the 
general solution is u — u reg + u S i ng , where 



Urcg = b n (z — 1)' 



and 



sing — 7 

n=0 



b n (z-l) n+1 ln(z-l) + d n (z-l) n ] 



(38) 



(39) 



The free parameters in the general solution are bo and dx- 
In the region where x < x c , z > 1. So by writing the equa- 
tion in terms of z we have found a singular solution whose 
logarithmic terms are real in the region where x < x c . 

When solv ing for corotating solutions we follow 
iBalbinskH l|l985f) and use the two different solutions in the 
regions where the logarithms have positive argument, and 
match at the corotation point. Denoting by ur the solution 
in the region where x > x c , and by ml the solution in the 
region where x < x c (and z > 1), we have 



UR 



UL 



Co + (Oo + ci)(x — X c ) 

+[oo + 01 (a; — x c )](x - x c ) In \x — x c \, 

d + (b +di)(z - 1) 

+[&o + 6i(2-l)](z-l)ln|2-l| 

d + (bo+di) {XC ~ X) 



(40) 



+ 



bo + fo 



(X c - X) 



111 



(41) 



In the limit as the solutions approach the corotation point, 
all of the terms tend to zero apart from the constant terms. 
For u to be continuous we must have 



co = do. 

Now consider the first derivatives. 



(42) 



du R 
dx 

dUL 

dx 



2a + ci + (3oi + 2c 2 )(a; — x c ) 

+ao In \x — x c \ + 2ai(x — x a ) In \x — x c \, 
z dun 
x c dz 



(43) 



Xc 
X 2 

-bo In 



2b + di + (36i + 2d 2 



(x c - X) 



261^^1 In 



(44) 



In the limit as the derivatives approach the corotation point 
all of the terms tend to zero apart from the constant terms 
and the plain logarithmic terms. We get 

du R 



lim 

x— >x~£ dX 

,. du L 
~~d — 



= 2ao + c\ + ao lim In la; — x c 



2bo + dx + bo lim In \x 



-bo lnx c 



(45) 



(46) 



Our insistence on the continuity of the function (equation 
1421 1 means that 



6 



-x c a . 



(47) 



Recall that ao and Co are related by equation 13411 . and so 
the coefficients of the singular In \x — x c \ term are equal 
at the corotation point. But in addition to this logarithmic 
term there is also a finite step in the first derivative. By per- 
mitting a step in the first derivative there will be a solution 
that meets the boundary conditions for any frequency within 
the range of the continuous spectrum. In fact, for each fre- 
quency there will be two solutions, one of each symmetry, 
with a particular step in the first derivative, that meet the 
boundary conditions (note that for a given frequency the 
step size for the solutions of opposite symmetry need not be 
the same). 

The step in the first derivative varies across the range of 
continuous spectrum frequencies and it can, for certain fre- 
quencies, be zero. If the step in the first derivative vanishes 
then 



2a + ci = (2&o + di — b inx c )- 

x c 

Using equation 1471 we can rewrite equation 1481 as 



di 



-x c (ci + ao lna;c). 



(48) 



(49) 



Removing the step in the first derivative means that 
there will be no delta function in the second derivative, 
merely a simple pole. Throughout the rest of this paper we 
will re fer to these solutions as zero-step solutions. iBal binski 
( 1985) suggested that such solutions were modes in the coro- 
tation region. Using the criterion of a vanishing step in the 
first derivative he was able to track the propagation of two 
such solutions into corotation as the amount of differential 
rotation was increased. These solutions subsequently merged 
to give birth to a dynamically unstable mode. Whilst it is 
not clear to us that these solutions are modes in the tradi- 
tional sense, they are clearly of interest if they can give rise 
to dynamical instabilities. They are also notable because at 
the corotation point the Wronskian of the two solutions be- 
ing matched (ur and ml) is zero. 

We have seen in this section that corotating solutions 
will in general be singular at the corotation point. At first 
sight this may seem unphysical. In solving the normal mode 
problem, however, it is implicitly understood that our ulti- 
mate goal is the evolution of initial data. The physical time- 
varying solution is obtained by convolving the solutions to 
the normal mode problem with the initial data and inte- 
grating over the frequency range (an inverse Laplace trans- 
form). The resulting integral solution will not be singular, 
and th e continuous spectru m will indeed have physical rele- 
vance JVan Kampenlll955l1 . In fact consideration of the ini- 
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tial value problem leads us to suspect that the zero-step 
solutions might exhibit behaviour that is physically distinct 
from the rest of the continuous spectrum. We are currently 
investigating this issue (Watts et al, in preparation). 



means that mQ, — a will contain a factor x 2 . As before we 
assume that 17' (0) 7^ and look for power series expansions 
of the form 



(53) 



7 SOLUTIONS ON THE BOUNDARY OF THE 
COROTATION REGION 

As discussed in Section [31 the r-modes are the only non- 
trivial solutions in the case of uniform rotation. An inter- 
esting question is whether these modes enter the corota- 
tion region as the degree of differential rotation is increased. 
In order to do this they must somehow develop corotation 
points. For the simple rotation laws we consider in this pa- 
per there are two possibilities. Either the corotation point 
appears first at the poles, or at the equator. In order for the 
mode to enter the corotation region, there must exist a solu- 
tion with a corotation point at the corresponding endpoint. 
We will consider each possibility in turn. 

7.1 Corotation point at the pole 

First consider (monotonic) rotation laws for which the angu- 
lar velocity is lowest at the pole. If modes are to enter coro- 
tation there must be mode solutions with corotation points 
at the pole. We can try to rule out such solutions by consid- 
ering the corresponding behaviour at the pole. This problem 
is easiest if we work with U, i.e consider equation 123H . 

If the corotation point is at the pole we can search for 
a power series expansion of the form 



The indicial equation yields 



U = a„ 



2\n + ~i 



(50) 



If fi'(l) 7^ then the indicial equation is 



1 4 2n'(i) K ' 

We now need to determine what values of 7 are per- 
mitted by the boundary conditions. It is clear from consid- 
eration of the physical variables that we require 5vg and 
Svp to be zero at the pole. From the continuity equation we 
can see that the second condition corresponds to requiring 
dU/dO = 0. Now U = Sv e sintf and 



dU _ 2^dU 

~dl ~ { )2 ~dx" 



(52) 



We therefore have two conditions to meet at the pole: (7 = 
and dU /dx less singular than (1 — a; 2 ) _1//2 . These conditions 
immediately rule out the possibility of 7 being imaginary. 
For real 7 the conditions exclude the possibility of 7 being 
negative, moreover we must have 7 > 1/2. Whether or not 
this condition can be met will depend on the rotation law 
and the degree of differential rotation. 



7.2 Corotation point at the equator 

Consider now rotation laws for which the rotation rate is 
lowest at the equator. If modes are to cross the corota- 
tion boundary there must be mode solutions with corotation 
points at the equator. Consider equation If the corota- 
tion point is at the equator the symmetry of the rotation law 



7 =2 



1± 1- 



8n'(o) 

Q"(0) 



(54) 



If 7 is real we have two possible power series, one for 
each value of 7 (which we call 7+ and 7-). 7+ > 1/2 and 
7_ < 1/2. The 7_ series has a singular derivative, which 
would correspond to a singular Sv v . If 7 were imaginary we 
would have 7 = 1/2 ±17/2, so once again the derivative 
would be singular. 

Suppose that we permit singular solutions but restrict 
the degree of singularity by requiring that the physical vari- 
ables Sve and 8v v be square integrable. This is a reasonable 
restriction as it enables us to define an energy measure. For 
this restriction to hold u' , and hence Sv v , must be less sin- 
gular than x~ x / 2 . The only permissible power series solution 
is therefore that with index 7+. 

If the maximum value of 7+ is less than 1 then tt'(0) 
is singular for even this solution. This would preclude the 
existence of even mode solutions with equatorial corotation 
points (for even modes we expect u'(0) = 0). In other words, 
it appears even r-modes cannot cross into corotation for cer- 
tain rotation laws. The development of equatorial corotation 
points in odd modes for these laws remains however a pos- 
sibility. 



8 APPLICATION TO SPECIFIC ROTATION 
LAWS 

The majority of the differential rotation laws examined in 
the literature were introduced because of their mathemat- 
ical simplicity rather than being motivated by the antici- 
pated physics. This is obviously fine as long as one is mainly 
interested in the properties of typical solutions, but it is 
clear that one ought to be somewhat cautious before draw- 
ing astrophysical conclusions from the relatively small set 
of rotation laws that have been considered to date. Ide- 
ally one would like to use a rotation law for a young neu- 
tron star that is grounded in physi cs, emerging for example 
from simulations of core collapse llZwerger fc Mulled Il997l : 
iDimmelmeier. Font fc Mullerll2002T) . At this stage, however, 
we are trying to deal with the technical challenges posed by 
the differential rotation problem. It therefore makes sense 
to consider simple rotation laws expressed in terms of a few 
parameters. Since we also wish to compare our results to 
previous work we have selected four rotation laws discussed 
in other papers. We have also considered a fifth law that 
exhibits two phenomena not seen for the other laws. 

The first two rotation laws, the j-constant and v- 
constant laws, h ave been used e xtensively in the literature 
(see for example lHachisul Ill986l0 . In particular, they were 
used in an examination of s table non-corotatine mode s of a 
differentially rotating shell jRezzolla fc Yoshidall200lll . The 
j-constant and v-constant laws are examples of rotation laws 
where the angular velocity is greatest at the pole. 
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The thi rd rotation law was discussed in a paper on solar 
oscillations (Wolff 1998). This is one of the few laws encoun- 
tered in the literature that is grounded in observations. Wolff 
used this law to examine oscillations of a differentially rotat- 
ing shell. This model caught our attention because Wolff's 
results suggested a large number of corotating solutions that, 
as differential rotation increased, combined with the more fa- 
miliar r-modes and disappeared, potentially giving rise to a 
dynamical instability. However, Wolff did not discuss either 
corotation points or dynamical instabilities. This law is also 
of interest as it is an example of a rotation law where the 
angular velocity is greatest at the equator. 

Our fourth rotation law is a simplified version of the 
Wolff law and i s a par ticul ar instance of a law discussed by 
both lBalbinskl l|l985h and lLuvtenl (|l990al) . This law is in- 
teresting because (as we will see) it exhibits both dynamical 
instabilities and regular solutions in corotation. 

The fifth rotation law has angular velocity greatest at 
the pole, but unlike the j-constant and v-constant laws, 
exhibits dynamical instabilities and modes with equatorial 
corotation points. 

All five laws obey the Solberg crit erion for dyn amical 
stability to axisymmetric perturbations jTassouil978Tl . This 
criterion, which generalises the Rayleigh condition for con- 
vective stability in an inviscid, incompressible fluid, states 
that we must have 

-^-(ft 2 m 4 ) > 0, (55) 
dzu 

where w, the cylindrical radius, is given by w = sin (9 = 
(f — x 2 ) 1 / 2 (we have set the radius of the shell equal to one). 

During violent events such as core collapse or neutron 
star neutron star merger it is conceivable that a star might 
achieve an intermediate rotation state that does not obey 
the Solberg criterion. The analysis conducted in the next 
section for the shell will only be valid for a Solberg unstable 
law if the timescale for growth of the Solberg instability is 
low compared to the timescales associated with the non- 
axisymmetric oscillations. 

8.1 Method 

In our analysis of each rotation law we took the following 
steps (in our view these would also be the natural steps in 
the analysis of other cases): 

(i) Check whether the law obeys the Solberg criterion 

(ii) Check whether the law has an inflexion point in the 
vorticity (fi' — 0). This will determine whether the law can, 
in principle, have dynamical instabilities or regular modes 
in corotation. 

(iii) Solve for real-frequency modes outside corotation. 
This is the easiest calculation to do, and it will indicate 
whether modes are approaching the corotation region. 

(iv) Solve for modes at the edge(s) of the continuous 
spectrum. This will provide information concerning whether 
modes cross into the corotation region. Furthermore, such 
solutions may (as we discuss later) indicate the onset of dy- 
namical instability. 

(v) Search for regular corotating modes 

(vi) Search for corotating modes with zero step in the first 
derivative at the co rotation p o int. T hese are solutions of the 
kind considered by Balbinskij lll985l) . 



(vii) Search for dynamically unstable modes developing 
at points where 

• Stable real-frequency modes breach the instability 
criterion 

• Stable real-frequency modes merge 

• Zero-step corotating solutions merge 

The notion that dynamical instabilities develop when 
modes merge is supported by empirical evidence from other 
instability calculations. The dynamical bar mode instability 
that limits the angular velocity of rotating stars sets in at 
points where two non-radial modes merge. The radial insta- 
bility that sets the maximum mass of a neutron star also 
sets in through the merger of two modes. Theoretical sup- 
port for this idea is given bv lSchuO (ll98Cf) . who notes that 
whenever instability sets in along a sequence the Hermitian 
eigenfrequency equation will have a double root. In addi- 
tion, we have already discussed how lBalbinskil jl985T) found 
dynamical instabilities developing at the point where two 
zero-step corotating solutions merged. 

The numerical integrations necessary to find the modes 
were carried out using a variable order variable step Adams 
method. A solution at one value of differential rotation was 
then used as an initial estimate for the solution at a slightly 
increased or decreased value. In this way we were able to 
track modes through the desired range of differential rota- 
tion. 

Modes without corotation points (real and complex fre- 
quencies) were located by picking a trial frequency, inte- 
grating from pole to equator and then iterating using the 
Newton- Raphson method until the equatorial boundary con- 
dition appropriate to the mode symmetry was met. A simi- 
lar routine was used to search for solutions on the corotation 
boundary, although in this case because the location of the 
corotation point was fixed there was a fixed relationship be- 
tween frequency and amount of differential rotation. 

Corotating solutions were found by splitting the domain 
of integration into two domains: x c < x < 1 and < x < x c . 
As the first domain has singularities at both ends we inte- 
grated away from both singularities and evaluated the Wron- 
skian at a matching point in the middle of the domain. On 
the second domain we integrated away from the singularity 
at x = x c and evaluated the function and its derivative at 
the equator. The desired solutions were those for which both 
the Wronskian on the first domain and either the function 
or its derivative at the equator (depending on the symmetry 
of the solution) vanished simultaneously. For both regular 
and general solutions in corotation the power series expan- 
sions of Section |H| were used to approximate the solution at 
x — x c ± S (where we used S ~ I0 -5 in the numerical calcu- 
lations). When searching for regular solutions in corotation 
for which Cl'(x c ) 7^ 0, we fix the arbitrary constant ao in 
the power series expansion and vary the frequency. For the 
general singular solutions the situation is a little more com- 
plicated. In the power series expansion for x — x c + S there 
are two free parameters, ao and ci . By fixing these and im- 
posing the matching conditions at x — x c we completely de- 
termine the two free parameters in the power series solution 
at x — x c — 8. In actuality the two desired functions become 
zero only for particular values of both the frequency and the 
ratio ci/ao (we use ci/ao rather than oo/ci as our variable 
because there are solutions where ci =0). We therefore fix 
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ao and vary both the frequency and this ratio, using a mod- 
ified version of the Powell hybrid method to locate zeroes of 
the two functions. Fortunately the frequencies and ratios of 
solutions turn out to be smooth functions of the parameters 
that determine the amount of differential rotation. 



8.2 Case 1: the j-constant rotation law 

The j-constant law for a shell takes the form 

n c A 2 



ft 



(56) 



1 + A 2 - x 2 ' 

The parameter A controls the amount of differential rota- 
tion , uniform rotation (ft = ft c ) being achieved in the limit 
as A — > oo. The rotation rate for this law is lowest at the 
equator. We get 

(57) 



^, _ 2ft c ,4 4 (l + A 2 + 3x 



(l + A 2 -x 2 ) 3 

Since ft' is always positive for real A there will be no dy- 
namically unstable modes or regular corotating modes for 
this rotation law. 

Figure shows the frequencies of a selection of modes 
and corotating zero-step solutions for m = 1 and m = 2, 
respectively. Note that for any value of A there is a singular 
solution at every frequency within the corotation band. In 
the Figures for this and the other rotation laws, however, 
we only show the zero-step solutions. 

The frequencies d erived for modes ou t side c orotation 
confirm the results of iRezzolla fc Yoshidal (j200lf> . As the 
amount of differential rotation increases, the frequencies of 
the r-modes decrease and approach corotation. The code 
crashes as a —> mft at the equator, crashing at larger values 
of A for modes that at uniform rotation have larger I. Modes 
that have very large I in uniform rotation are not shown, but 
follow the general trend illustrated in the figures. Whether 
modes continue to exist very close to the corotation bound- 
ary as the degree of differential rotation is increased still 
further is difficult to say. We have however tracked modes 
very close to the corotation boundary and our results suggest 
that modes outside corotation cease to exist below a certain 
value of A. Disappearance of modes as a system parameter 
changes is not unexpected. Several examples of eigenvalue s 
vanishing when parameters change are given in lKatol 11995). 

The question of whether modes can enter the corota- 
tion region by developing equatorial corotation points was 
discussed in section For this law the maximum value of 
7+ is 1. This means that, analytically, we were able to rule 
out any even solutions with equatorial corotation points for 
any value of A, and any odd solutions for A > 0.354. This 
was confirmed by our numerical results. For A < 0.354 we 
were not able to rule out analytically the possibility of odd 
r-modes crossing the boundary. Our numerical results sug- 
gest however that there are no odd r-mode solutions existing 
below A w 0.36. For this rotation law, therefore, we find no 
evidence that r-modes cross into corotation. 



1 In the limit A — > the j-constant law tends to zero rotation. We 
do not discuss this limit further in this paper and hence assume 
A > 0. 



8.3 Case 2: The v-constant rotation law 

Another law that has been examined in the literature is 
the v-constant rotation law. There are however two differ- 
ent laws_with_this na me in the literature . The first, used 
bv lHachisul Jl986ft and lKarino fc Eriguchl J2002I) . takes the 
form 



ft 



Q C A 



(58) 



(A 2 + i-x 2 y/ 2 ' 

The second, used by lEriguchi fc MiilleJ 

<1985Il iKarino Yoshida fc Eriguchil (|200Jl . and 
IRezzolla fc Yos hida (2001), takes the form 



ft = 



il c A 



a + (l- x 2 y/ 2 ' 



(59) 



In both cases uniform rotation (ft = ft c ) is reached in 
the limit A — > oo. For the first law it is not possible to have 
ft' = 0, so we can rule out the possibility of dynamically 
unstable modes or stable modes within corotation. For the 
second, ft' can be zero for A > 1 or A < 1/2, so we cannot a 
priori rule out the possibility of unstable modes. No instabil- 
ities were found, however, by our numerical analysis. In fact 
the modes and continuous spectrum for the two v-constant 
laws behave qualitatively identically to the modes and con- 
tinuous spectrum of the j-constant law. For this reason we 
have not included figures of the results. 



8.4 Case 3: Wolff's solar rotation law 

Wolff's solar rotation law fWblfl1ll998h takes the form 



ft = 2tt(454.8 - 60A/3x 2 - 71.4/fe 4 



(60) 



The parameter /3 measures the amount of differential rota- 
tion, uniform rotation corresponding to f3 = 0. Wolff consid- 
ers the range < (3 < 1 but we will consider an upper limit 
of /3 = 2 as there are interesting phenomena in this range. 
The rotation rate for this law is lowest at the pole. We get 



ft' = 2tt(909.6 + 120.8/3 + 132/3a; 2 - 2142/3a; 4 



(61) 



In this case ft' has a zero on the domain < x < 1 for f3 > 
0.4815. Both dynamical instabilities and regular solutions in 
corotation are therefore possible. 

Figure |21 shows the frequencies of a selection of modes 
and zero-step corotating solutions for m = 1 and m — 2. As 
the amount of differential rotation increases, the frequency 
of the each r-mode shifts and in due course all modes en- 
counter the corotation boundary. Modes of low I appear to 
approach the boundary directly, modes of higher I seem to 
approach tangentially. Searching for mode solutions on the 
corotation boundary, as described in section 17.11 we find 
valid solutions at the points where the I = (2,3),m = 1 
and I = (2, 3, 4), m — 2 modes encounter the boundary, sug- 
gesting that these modes can develop corotation points 2 . As 
m increases so does the number of valid solutions on the 
corotation boundary. For m = 3 and m — 4 there are valid 



2 Figure[5]appears to suggest that the I = 4, m = 1 mode crosses 
the boundary. There is however no physically valid solution on 
the boundary. Resolving this region in more detail it is possible 
to see the I = 4 mode and the corresponding zero-step corotating 
solution approaching the boundary tangentially but not crossing. 
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Figure 1. Results for rotation law 1: A determines the degree of differential rotation, rotation being uniform in the limit A — ► oo. The 
plot on the left shows m = 1, that on the right m = 2. Solid lines: real-frequency modes outside corotation. Dashed line: The lower 
boundary of the corotation region - solutions to the right of this line are in corotation. Dotted lines: real-frequency solutions with zero 
step in the first derivative at the corotation point. The values of I given refer to the uniform rotation limit A — > oo, in which, as discussed 
in section |21 there is one r-mode solution for each value of I. 
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Figure 2. Results for rotation law 3: f3 determines the degree of differential rotation, rotation being uniform when /3 = 0. The plot on 
the left shows m = 1, that on the right m = 2. Solid lines: real-frequency modes outside corotation. Dashed line: The lower boundary of 
the corotation region - solutions to the right of this line are corotating. Dash-dot lines: Real part of frequency for dynamically unstable 
modes. Dotted lines: real-frequency solutions with zero step in the first derivative at the corotation point. For dynamical instabilities to 
exist, we require /3 > 0.4815. The values of I given refer to the uniform rotation limit /3 = 0, in which, as discussed in section[3] there is 
one r-mode solution for each value of I. 



solutions at the point where the four lowest I modes en- 
counter the corotation boundary. Modes of larger I, which 
appear to approach tangentially, do not seem to cross the 
boundary. 

A numerical search for modes that are regular in coro- 
tation but for which the corotation point does not coincide 
with the zero of Q' proved fruitless. For the range of f3 exam- 
ined there were no valid solutions on the domain x c < x < 1 
(see the discussion in Section fS.H . 

What about the possibility of regular solutions, with 
corotation points at the point where f2' = 0? Such solutions 
could only exist at one frequency for a given f3, this frequency 
being a smooth function of j3. If we represent the solution 
as a sum over spherical harmonics, U = CpY™, it is 
possible to derive a recursion relation for the coefficients C; m . 
A numerical search for these modes yields nothing, however, 
and a numerical examination of the coefficients C ; m shows 
why. The coefficients are divergent, with CJ™ — > oo as I — » oo. 



There are therefore no regular modes in corotation for this 
law. 

Continuous spectrum solutions with zero-step in the 
first derivative are present for this rotation law. Such solu- 
tions are generated at the point where real-frequency modes 
appear to enter corotation (i.e. where there are valid solu- 
tions on the corotation boundary). There are also many such 
solutions present at low differential rotation which, in con- 
trast, seem to approach both edges of the corotation region 
tangentially. 

Dynamical instabilities are found to develop at the 
point where modes encounter the corotation boundary at 
values of f3 for which Cl' has a zero on the shell (in accordance 
with the instability criterion derived in Section [KJ. Modes 
that encounter the corotation boundary below the threshold 
value of f3 do not exhibit instabilities. There are five unstable 
modes, corresponding to / = (2, 3),m = 1;Z = (2, 3)m = 2 
and I = 3, m = 3. We find no instabilities for modes of 
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higher I or m. Frequencies and growth times for these un- 
stable modes are shown in figure [3] The behaviour of the 
I = m — 3 mode is particularly interesting as it indicates 
that the mode becomes stable again for high differential ro- 
tation. The I = 3, m — 2 mode exhibits similar behaviour at 
higher j3. 

We can see from equations 1251 and 1611 why we should 
not expect to find any unstable modes at high (3. As f3 in- 
creases, the negative part of ti' will come to dominate the 
integrand in equation 12511 . Beyond a certain threshold value 
it will no longer be possible to satisfy the integral equation. 
The precise value of the threshold will depend on U but all 
modes will reach this point for some f3. 

The results reported in this paper agree with those re- 
ported by Wolff for real-frequency modes outside corotation. 
However, they differ dramatically for solutions with corota- 
tion points. In addition, we find dynamical instabilities not 
mentioned by Wolff. The discrepancy between the results 
likely originates with Wolff's perturbation method and with 
his choice of basis eigenfunctions. He uses not only the stan- 
dard r-mode eigenfunctions discussed in section |3] but a set 
of toroidal geostrophic eigenfunctions of single l,m. But the 
geostrophic solutions form a degenerate set in uniform rota- 
tion and cannot be classed as eigenfunctions. As such they 
cannot be used as a basis for perturbation. 

We can compare the results from our investigation with 
measurements of the deg ree of differential rotation at var- 
ious depths in the Sun dHowe et a The degree of 
differential rotation is highest in the surface layers and the 
convective envelope, where values of Q(0)/2n = 450 nHz 
and f2(cos(7r/6))/27r = 370 nHz are reported. This would 
correspond to having (3 ~ 0.9. At this level in our model we 
see dynamically unstable modes, the fastest having a growth 
time of just a few years. 

Our model is very simplistic. It does not, for exam- 
ple, include viscosity, which would be expected to shift the 
points of marginal stability to a higher degree of differen- 
tial rotation. However, in the outer layers of the Sun the 
degree of differential rotation is almost independent of the 
radius, suggesting that a shell model may not be a bad ap- 
proximation. In this case a dynamical instability of the type 
found in this model could be involved in limiting the de- 
gree of differential rotation observed in the Sun. See also 
iDzhalilov. Staude fc Oraevskvlll200 2'l for a discussion of var- 
ious Solar properties that might be explained by the action 
of unstable r-modes with growth times of ~ 10 — 1000 years. 
Whether dynamical instabilities of the type found in this pa- 
per could play a similar role is an open question. 



8.5 Case 4: A simple quadratic rotation law 

The fourth rotation law we have considered takes the form 



n = fi c (i - fix 2 



(62) 



Here the parameter f3 measures the degree of differential ro- 
tation and we have examined the range < j3 < 1. Uniform 
rotation (fi = Q c ) is recovered if f3 = 0. As with Wolff's law 
the rotation rate for this law is greatest at the equator. We 
get 



From this we see that Cl' has a zero on the domain if f3 > 0.2. 
Both dynamical instabilities and regular solutions in corota- 
tion are therefore possible. Figure^Jshows the frequencies of 
a selection of modes and zero-step corotating solutions for 
m — 1 and m = 2. 

Just as for Wolff's law, modes of low I approach the 
corotation boundary directly as differential rotation is in- 
creased. Modes of higher I tend to approach tangentially. 
And in accordance with this there are valid solutions on 
the corotation boundary for the I = (2,3),m = I and 
I = (2, 3), m = 2 modes. 

The I = 3 mode is a little different however, as it passes 
straight into corotation as a regular mode. To see why this 
happens consider the following argument. For this rotation 
law we can write 



mf) — g 
and 

$V = -12/?fic 



-m(3Q, c 



2 (mflc — a) 



mQ, c (3 



1 + /3 



6 i 



(64) 



(65) 



As discussed in sect ion lo"Tl the zeroes of these two functions 
may coincide, making the governing equation regular and 
giving rise to regular solutions in corotation. This will occur 
if 



mfic(5-/3) 
" = 6 ' 

In this case the governing equation 1231 becomes 

m 2 U " 



(I - x 2 )U" - 2xU' 



l-x 2 



+ 6(7 = 0. 



(66) 



(67) 



This equation is simply Legendre's equation for 1 = 3. Thus 
there is a pure 1 = 3 mode that corotates but is completely 
regular. Just as for Wolff's law, however, there are no regular 
corotating solutions of the type given in equation 1281 . 

Only the I = 2, m = 1 and I = 2, m = 2 modes become 
dynamically unstable at the point where they encounter the 
corotation boundary. These are the only modes that hit the 
boundary after passing the threshold value of f3 for which JY 
first has a zero on the shell. Frequencies and growth times 
for both modes are shown in figure |S] 



8.6 Case 5: A simple trigonometric rotation law 

The final rotation law we have considered is 



fi = Q c (1 — f3 cos x) 



(68) 



O' = 2fi c (l + {3-6f3x 2 



(63) 



The parameter (3 measures the degree of differential ro- 
tation and we have examined the range < (3 < 1. Uniform 
rotation (Q = f2 c ) is recovered if (3 = 0. The rotation rate for 
this law is greatest at the pole, as it is for the j-constant and 
v-constant laws. What is different about this law, however, 
is that Cl' = for j3 > 2/3, which could permit dynamical 
instabilities. Figure |S| shows mode frequencies for m = 2 for 
a reduced range of the differential rotation parameter (3. The 
figure shows interesting phenomena that we do not see for 
any of the other rotation laws. Modes crossing the corota- 
tion boundary for this law would have to develop equatorial 
corotation points, as discussed in section For this law, 
the maximum value of 7+ is 2. Modes of odd symmetry can 



Differentially rotating spherical shells: Normal modes 11 





a/£2c a/flc 



Figure 4. Results for rotation law 4: f) determines the degree of differential rotation, rotation being uniform when /3 = 0. The plot on 
the left shows m = 1, that on the right m = 2. Solid lines: real-frequency modes outside corotation and regular modes inside corotation. 
Dashed line: The lower boundary of the corotation region - solutions to the right this line are corotating. Dash-dot lines: Real part of 
frequency for dynamically unstable modes. Dotted lines: real-frequency solutions with zero step in the first derivative at the corotation 
point. For dynamical instabilities to exist we require (3 > 0.2. The values of / given refer to the uniform rotation limit /3 = 0, in which, 
as discussed in section [3] there is one r-mode solution for each value of I. 



in principle cross the boundary if (5 > 16/25, modes of even 
symmetry if j3 > 2/3. 

Let us start by considering the I — 3 mode. As the 
amount of differential rotation is increased, the frequency 
approaches the corotation boundary. The mode crosses the 
corotation boundary above the instability threshold value of 
/?, and above this point continues to exist as a dynamically 



unstable mode. An odd zero-step solution also develops at 
this point, but this is not shown in Figure |S] as the fre- 
quency is very close to the real part of the frequency of the 
dynamically unstable mode. 

The situation for the even modes is slightly different. 
Although the modes could in principle cross the corotation 
boundary, we find that the I = 2, m — 2 and I = 4, m = 
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Differential rotation (beta) Differential rotation (beta) 

Figure 5. Dynamically unstable modes for law 4. The left frame shows the real part of the frequency in dimensionless units (Re<r/n c ). 
The right frame shows the growth time in dimensionless units (£7 C /Imcr). 



2 modes merge outside corotation, becoming dynamically 
unstable. Such a route to instability is not seen for the other 
rotation laws examined. With reference to Figure HJ note 
that we were unable to track the I = 4 mode for 0.644 < (3 < 
0.678 because the frequency was too close to the corotation 
boundary. We are however confident that it is this mode 
that is merging with the I — 2 mode, for three reasons. 
Firstly, the merging mode has even symmetry (as we would 
expect). Secondly, all even modes with I > 4 approach the 
corotation boundary at lower values of (3. The I — 4 mode is 
the only other even mode apart from I = 2 that appears to 
persist to high j3. Thirdly, the merging mode is not emerging 
from corotation, as the only valid solution on the corotation 
boundary for m = 2 occurs at the point where the I — 3 
mode crosses. The most reasonable conclusion is that the 
I — 4 mode continues to exist just outside corotation for 
0.644 < (3 < 0.678 before moving away from the boundary 
to merge with the / = 2 mode at higher j3. 



9 SUMMARY 

The oscillations of neutron stars are thought to be a promis- 
ing source of gravitational waves. In order to develop accu- 
rate models of gravitational wave emission we need to de- 
velop a full understanding of the dynamics of these stars. 
One factor that has been neglected in many previous stud- 
ies is differential rotation. Neutron stars are thought to be 
born with differential rotation, and there are several plausi- 
ble mechanisms that could sustain differential rotation be- 
yond birth. Differential rotatio n introduces the poss ibility 
of dynamical shear instabilities iPapaloizou fc Pringlelll984l 
1985). In addition, the perturbation equations of differen- 
tially rotating inviscid systems are complicated by the pres- 
ence of singularities at points where the pattern speed of a 
mode equals the local angular velocity of the background 
flow. These singularities mean that differentially rotating 
inviscid systems possess a continuous spectrum in addition 
to discrete modes. To begin to understand this problem we 
have in this paper analysed the normal mode problem for a 
toy model, a differentially rotating shell. Although the model 



is simplistic, it exhibits many of the key features of the more 
complex three-dimensional problem. 

In the uniform rotation limit the shell problem admits 
the standard r-modes as solutions. As the degree of differen- 
tial rotation increases the frequencies of these modes grad- 
ually approach the lower edge of the expanding corotation 
band (m,fl m i n < a < mfl ma x)- Modes that cross the corota- 
tion boundary are found to be rare. Indeed, we can rule out 
completely the crossing of even modes for certain rotation 
laws. For other rotation laws and symmetries we find cross- 
ings only for modes that have low / in the uniform rotation 
limit. All other modes approach the boundary tangentially 
and appear to cease to exist when the degree of differential 
rotation exceeds a threshold value. Modes that have higher I 
in the uniform rotation limit cease to exist at lower degrees 
of differential rotation. 

For frequencies in the corotation band the shell prob- 
lem admits a continuous spectrum of corotating solutions 
with discontinuous first derivatives. Although such solutions 
may appear at first sight unphysical, they will have physical 
significance when we consider the initial value problem. In 
general, corotating solutions have both a logarithmic singu- 
larity and a finite step in the first derivative at the coro- 
tation point. For certain special frequencies, however, the 
step in the first derivative vanishes. In differentially rotating 
cylinders these zero-ste p solutions can m erge to give rise to 
dynamical instabilities jBalbinskil ll 985 'l . but we observe no 
such behaviour in our problem for the rotation laws exam- 
ined. We cannot rule out the possibility, however, for other 
rotation laws. 

The character of the modes that cross into corota- 
tion varies. At all points where modes cross the corotation 
boundary we observe the development of zero-step corotat- 
ing solutions (except in one unique case, discussed below). 
At some points we also observe the development of dynami- 
cally unstable modes. A necessary (but not sufficient) condi- 
tion for dynamical instability is that Cl' be zero somewhere 
on the shell. For the rotation laws examined this condition is 
met only when the degree of differential rotation exceeds a 
threshold value. All modes crossing the corotation boundary 
above this threshold become dynamically unstable. We also 
observe, for one rotation law, the development of a dynam- 
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Figure 6. Results for rotation law 5 for m = 2: f3 determines the degree of differential rotation, rotation being uniform when /3 = 0. 
Solid lines: real-frequency modes outside corotation. Dashed line: The lower boundary of the corotation region - solutions to the right of 
this line are corotating. Dash-dot lines: Real part of frequency for dynamically unstable modes. Dotted lines: Real frequency solutions 
with zero step in the first derivative at the corotation point. For dynamical instabilities to exist we require /3 > 2/3. The values of I given 
refer to the uniform rotation limit j3 = 0, in which, as discussed in section [3] there is one r-mode solution for each value of I. 



ical instability via the merger of two stable modes outside 
corotation. Only in one exceptional case for the simplest ro- 
tation law did we observe a mode crossing the corotation 
boundary and continuing to exist as a regular stable mode. 
Such regular corotating modes are unlikely to exist for more 
realistic rotation laws. 

Before moving on to the three-dimensional problem, we 
intend to investigate two related problems for the simple 
shell model. The first factor of interest is the effect of vis- 
cosity, an essential component of any realistic model. If we 
consider the Navier-Stokes equations instead of the Euler 
equations the vorticity equation for the stream function be- 
comes 



(a-mQ)V 2 e U + 



mdeCl 
sin 6 



U 



Ve(y 2 e U) + 



(sin 9 cos 9dgU — (cos 9 + m )Uj 



0. 



(69) 



where v is the coefficient of shear viscosity. Equation 1691 1 
is not singular, as the coefficient of the leading derivative is 
never zero. So in the viscous problem there is no continuous 
spectrum, simply sets of discrete modes. In the zero viscosity 
limit, however, the solutions to the viscous problem should 
reduce to the modes and continuous spectrum of the invis- 
cid problem. At first sight it is difficult to see how discrete 
well-behaved modes could reduce to the singular continuous 
spectrum solut ions that we derived for the inviscid problem 
m Section El iBovdl <llQ8 lT) resolves this difficulty by not- 
ing that any physical continuous spectrum perturbation is 
made up of an integral over frequency of the singular invis- 
cid solutions. When solving the viscous problem, the modes 
are discrete for even the tiniest viscosity and no integration 
over frequency is required. Solution of the viscous problem 
should therefore shed light on the physical manifestation of 
the continuous spectrum. 

We are also pursuing the initial value problem. It is 
clear from our analysis of the normal mode problem that 
the time dependence of the continuous spectrum will only 
be understood by solving the initial value problem. In fact 



the singular continuous spectrum solutions found in this pa- 
per only acquire physical meaning when we consider the 
evolution of initial data. We are particularly interested to 
see whether the zero-step solutions appear physically dis- 
tinct in time evolutions. That the zero-step solutions might 
exhibit different time-dependence is suggested by the fact 
that these solutions have zero Wronskian, which will con- 
tribute an additional pole to the integrand of the Laplace 
Transform inversion integral (Watts et al, in preparation). 
We also hope that an analysis of the initial value problem 
will enable us to develop necessary and sufficient conditions 
for instability. 
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